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Abstract. Pulsations in subdwarf B stars are attributed to radiative levitation of iron- 
group elements in the stellar envelope. Until now, only iron diffusion is accounted for 
in stellar models used for sdB seismology. However, nickel has also been suggested 
as a contributor to the opacity bump that drives the pulsation modes. Stellar models 
including time-dependent atomic diffusion, as we compute here, are needed to evaluate 
the importance of different iron-group elements for mode driving. We perform detailed 
calculations of radiative accelerations of H, He, C, N, O, Ne, Mg, Fe and Ni and include 
these in Burgers' diffusion equations. We compute the evolution and non-adiabatic pul- 
sations of a typical subdwarf B star We show that, despite its lower initial abundance, 
nickel accumulates to comparable mass fractions as iron in the sdB envelope. For ac- 
curate determination of pulsation frequencies and mode instability, it is essential that 
diffusion of both metals are included in stellar models. The role of other iron-group 
elements remain to be evaluated. 



1. Introduction 



Subdwarf B (sdB) stars have been identified as core He burning stars with an unusual 
thin H-envelope (< 0.02 M©), also called Extreme Horizontal Branch stars because of 
their location in the Hertzsprung-Russel diagram ( Heber 1986) . Two subgroups of sdB 
stars are pulsating with variable class names V361 Hya and VI 093 Her. The discovery 
of short-period (100^00 s) pulsations in V361 Hya stars by 'Kilkenn y et al.] ( \991) 
coincided with the prediction of unstable pressure(/7)-mode pulsations in these stars by 
Charpinet et al. ( 1996| l. The excitation was attributed to the opacity mechanism enabled 
by Fe accumulating diffusively in the stellar envelope ( Charpinet et al.|[T997 1. A few 



years later, Green et al. 



(20031 discovered long -period pulsations (30-120 min) in the 



Fontaine et al.| ( |20()3 1 showed that the same opacity mechanism 



cooler V1093 Her stars. 

can also excite long-period gravity(g)-modes. However, these authors found a too cool 
theoretical blue-edge of the g-mode's instability strip compared to observations. The 



problem can be partly solved by using OP opacities ( [Badnell et al. 2005 1 instead of 
op al's (jlglesias & Rogers^ '1996) and by assuming that Ni accumulates as well as 
Fe (Jeffery & Saio 20061. Furthermore, Hu et al. ( [2009 , ) showed that by including 
H-He diffusion the theoretical blue-edge is shifted even closer to the observed value 
of 30 kK. Still, all aforementioned models fail to produce unstable g-modes at low 
spherical degree {€ < 3) at the observed temperatures, while these modes are the most 
likely to be observed owing to geometric cancellation of higher degree modes. It is 
clear that stellar models with time-dependent diffusion, such as we produce here, are 
needed in sdB seismology. 
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2. Method 



We use a version of the stellar evolution code stars ( Eggleton|1971 1 that is adapted for 
asterose i smolo gy ( |Hu et al.||2008[ ). The code uses nuclear reaction rates from [Angulo 
&etal.|(|l999'), except for the~'^(p, 7)^^0 rate for which the recommended value by 



Heryyig et al. (2006 ) and [Formicola & LUNA Collaboration (2002) is used. Neutrino 
loss rates are according to Itoh et al. ( 1989 1992 ). Convection is treated with a standard 
mixing-length prescription (Bohm-Vitense 1958) with a mixing length to pressure scale 
height ratio of 1/Hp = 2.0. Convective mixing is treated as a diffusive process in the 
framework of mixing length theory. The Rosseland mean opacity a-r and radiative 
accelerations ^rad are computed in-line during the evolution to account for composition 
changes consistently. For this we use atomic data and codes from the Opacity Project 
(OP, [Seaton et al.| 1994| and |Badnell et al.| | 2005|l. For high ternperat ures outside the 
OP range {T > 10 K), we use OPAL tables ( Iglesias & Rogers]] 1996 1 combined with 
conductive opacities ( Cassisi et al.|2007" l 



The starting model for each of our calculations is a ZAEHB model with total 
mass M* - 0.46, envelope mass Menv - 3.5 x \0~'^Mq and metallicity Z = 0.02 with 



Grevesse & Noels] ( [1993] ) metal mixture. We follow the stellar evolution from the start 
to the end of of core He burning. We compute two series of simulations with atomic 
diffusion of H, He, C, N, O, Ne, Mg and 

(a) Fe; 

(b) FeandNi. 



The diffusion velocities are computed by solving the Burgers ( 1969 1 equations with 



a version of Thoul et al. ( 1994) routine that is adapted to use Paquette et al. (19861 
resistance coefficients, mean ionic charges and radiative accelerations. 

We use the linear non-adiabatic pulsation code mad (Dupret' 2001") to compute 
the seismic properties, in particular pulsation frequencies and damping rates, for the 
stellar evolutionary models. We restrict ourselves to modes of spherical degree € <2, 
because higher £ values are geometrically disfavoured for observation. We search for 
eigenfrequencies in the range of 0.1 < w < 20 where u is the dimensionless angular 
eigenfrequency, cd - InfTdyn, f is the frequency and Tdyn is the dynamical time-scale. 

The method summarized here is described in more detail by |Hu et al (201 1 1. 



3. Results 



We show the effects of atomic diffusion on the stellar evolutionary track in Sect. 3. 1 In 



Sect. 3.2 the changes in abundances and opacities as result of diffusion are discussed. 



The effects on mode driving are studied in Sect. 3.3 



3.1. Evolutionary tracks 

In Fig. [T] we plot the evolutionary tracks of the two simulations in a Tgff - log g diagram. 
We see that including Ni diffusion causes a shift in the track compared to only Fe 
diffusion. The star becomes slightly cooler and more extended due to an increase in 
opacity (see Fig.[4]l. The pulsation frequencies, in particular those of ;?-modes, are very 
sensitive to log g and it is therefore important to model diffusion accurately. 
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Figure 1 . Evolutionary tracks in the T^s - log g plane for sdB stars with atomic 
diffusion. Simulation (a) is with diffusion of H, He, C, N, O, Ne, Mg and Fe; simu- 
lation (b) is with diffusion of H, He, C, N, O, Ne, Mg, Fe and Ni. The crosses mark 
the part of the evolution for which the pulsations are shown in Fig.|5] 



3.2. Abundances and opacities 

In Fig. |2] and Fig. |3] the time evolution of the interior abundances of Fe and Ni can 
be followed for simulation (a) and (b), respectively. We see that, despite its lower ini- 
tial abundance, Ni accumulates to comparable mass fractions as Fe in the sdB envelope. 
Further, whether or not Ni is allowed to diffuse has an effect on the way Fe accumulates. 
This is because the chemical composition of the mixture influences the radiative accel- 
eration of each element. After all, the radiative acceleration depends on the Rosseland 



mean opacity and the monochromatic opacity of the mixture (see e.g. |Seaton|1997^ . 

The Rosseland mean opacity profiles are shown in Fig. [4] It can be seen that 
Ni accumulation causes the opacity bump to occur at slightly higher temperatures as 
akeady pointed out by Jeffery & Saio ( |2006[ l. 



3.3. Stability analysis 

In Fig. [5] we show the pulsation periods as a function of T^^. We plot pulsations for only 
a part of the evolution (indicated in Fig. [T]l so that the graphs do not become overclut- 
tered. We see that Fe diffusion alone cannot efficiently drive ^-modes. By including Ni 
diffusion, many g-modes become unstable but the driving of /7-modes becomes slightly 
less efficient. This is related to the location of the iron-group opacity bump as shown in 
Fig.|4]and discussed by |JefFery & Saio (2006 1. It is worth emphasising that we are plot- 



ting modes of degree £ <3. These are the first sdB models that show unstable ^-modes 
at such high effective temperatures (Jeff ~ 29 kK). In fact, the blue-edge problem of the 
V1093 Her instability strip can finally be solved with our models (Hu et al. in prep.). 
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Figure 2. Interior abundance profiles of Fe and Ni at different EHB ages for sim- 
ulation (a). 
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Figure 3. Interior abundance profiles of Fe and Ni at different EHB ages for sim- 
ulation (b). Lines are the same as in Fig.|2] 
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Figure 4. Rosseland mean opacity profiles as a function of temperature at different 
EHB ages. 
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Figure 5. Pulsation periods as a function of effective temperature. Red dots are 
unstable modes and green crosses are stable. The left panel is for simulation (a) and 
the right panel for simulation (b). 
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4. Conclusions 

We have shown the importance of Ni diffusion, in addition to Fe diffusion, for the study 
of sdB pulsations. Not only is there a large effect on the mode driving, the stellar 
evolutionary track is also affected by Ni diffusion. This is expressed in shifts in T^ff 
and log g and hence the pulsation frequencies. An accurate prediction of the theoretical 
pulsation frequencies is essential now that space data from CoRot and Kepler provide 
high-precision observed frequencies. It may be necessary to treat diffusion of other 
(iron-group) elements as well, in particular Cr and Mn, and we shall explore this in a 
forthcoming paper. 
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